All-magnonic spin-transfer torque and domain wall propagation 
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The spin-wave transportation through a transverse magnetic domain wall (DW) in a magnetic 
nanowire is studied. It is found that the spin wave passes through a DW without reflection. A 
magnon, the quantum of the spin wave, carries opposite spins on the two sides of the DW. As a 
result, there is a spin angular momentum transfer from the propagating magnons to the DW. This 
magnonic spin-transfer torque can efficiently drive a DW to propagate in the opposite direction to 
that of the spin wave. 
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Magnetic domain wall (DW) propagation along 
nanowires has attracted much attention in recent years 
[D-Q because of its fundamental interest and potential 
applications [HQ- So far, a spin-polarized electric cur- 
rent and/or magnetic field including a microwave 3 are 
the two known control parameters for manipulating DW 
propagation along nanowires: A magnetic DW propa- 
gates along a wire under a static magnetic field because 
of energy dissipation Q while a DW moves under an elec- 
tric current because of spin- transfer torque (STT) UGH. 
In terms of spintronic applications based on an electron 
spin current STT, the Joule heating due to the excessive 
high critical current density @, | is a bottleneck. Thus, 
it should be very interesting and important both academ- 
ically and technologically if one can find other effective 
control methods and principles for DW manipulation in 
magnetic nanowires. 

Both electrons and magnons, quanta of spin waves, 
carry spins. A magnon is a spin-1 object with an angular 



momentum of h 
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Similar to the STT from electrons 
to magnetization, a STT from magnons to magnetization 
should in principle exist. Indeed, polarized electric cur- 
rent generated by heat induced magnons in spin valves 
was predicted theoretically [l2|, [l3| and was confirmed 
experimentally 1A\ . The interaction between spin waves 
and DW had also been investigated quantum mechan- 
ically 15( and classically 16H18|. The time-dependent 
Schrodinger equation was used [15( to show that a DW 
is stable when it interacts with a spin wave, and the spin 
wave is reflected by the DW which is different from our 
finding below. Of course, they did not study magnonic 
STT. However, the phase change of spin waves after pass- 
ing through a DW was predicted. In terms of STT di- 
rectly from magnons, the question is how one can facil- 
itate a spin exchange between magnons and magnetiza- 
tion. In this Letter, we show that the spin wave inside a 
DW satisfies a Schrodinger equation with a reflectionless 
potential well. A magnon changes its spin by 2H (the 
magnon spin flips from — H to H) after passing through 
the DW, as shown in Fig. 1. This angular momentum is 
absorbed by the DW, resulting in propagation of the DW 
in the opposite direction to that of spin-wave propaga- 



tion. The validity of these findings is verified by solving 
the Landau-Lifshitz-Gilbcrt (LLG) equation numerically 
in a one-dimensional nanowire with material parameters 
of ferrimagnet yttrium iron garnet (YIG). The frequency 
and the field dependences of the DW propagation speed 
are also studied. 

Consider a head-to-head DW in a magnetic nanowire 
whose easy axis defined as the z axis is along the wire as 
shown in Fig. 1, the magnetization dynamics is described 
by the LLG equation [J] , 



dm 

— - = m x h cff 
at 



dm 

amx— , (1) 



where m is the unit direction of local magnetization 
M = mM s with a saturation magnetization M s , a is the 
phenomenological Gilbert damping constant, and h ff is 
the effective magnetic field consisting of anisotropy and 
exchange fields in the unit of M s . t is normalized by 
(7M S ) _1 and 7 is the gyromagnetic ratio. For simplicity, 
we consider a uniaxial wire with h e g = Km z z + A^^ 
where m z is the z component of m, and K and A are the 
anisotropy and exchange coefficients, respectively. In the 
spherical coordinates of polar angle 9 and azimuthal an- 
gle <fi, m = (sin 9 cos </>, sin 9 s'mcj), cos 9). For a static DW, 
m = mo is given by the Walker profile tan ^ = exp (^-) 
and lies in a fixed plane, say the y-z plane (4>q = 7r/2), 
where A = ^/A/K is the DW width Q. 

To derive the equation of motion for the spin wave, a 
small fluctuation of m around mo is expressed in terms 
of unit directions e r , (■■„, and e, ; , defined by mo- 



m = e. 



+ [m e (2) e g + m^ (z) e^] e" 



(2) 



where u is the spin- wave frequency, mg and are small, 
m 2 e + m 2 < 1. Substituting Eq. ([2} into Eq. ^ and 



neglecting the higher-order terms, such as mg, mgm' g , 
mgm^, etc., (' denotes the derivative in z), we obtain, in 
the absence of the damping, 

— iuimg = Am'^ + K (2 sin 2 9 — l) m^, (3) 
iujm,p = Am'g + K (2 sin 2 #0 — l) fne- (4) 
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FIG. 1: (Color online) Illustration of a transverse DW struc- 
ture whose m is denoted by the (blue) arrows. The spin wave 
(magnon) is a small amplitude precession of m (represented 
by the red cones) around the static DW. A linearly polar- 
ized microwave h (t) is applied in a small region on the left 
side of the DW so that the generated spin wave propagates 
through the DW from the left. A is the DW width. Inset: The 
magnons (wavy lines with arrows indicating the propagating 
directions) pass through the DW (represented by the rectan- 
gular block) from the left to the right without any reflection. 
The magnon spin (solid circle with an arrow) is —h on the left 
side of the DW and h on the right side. The magnonic STT 
drives the DW propagating to the opposite direction (green 
arrow) of the spin wave, with the velocity Vdw- 



Defining tp = mg — irrij,, ([3]) and Q can be recasted as 



«M0 = 



-, and q 2 



(5) 



with £ 

equation with propagating waves [20|, |2l| , 



1. This is a Schrodingcr 



v(0 = p 



tanh£ 



-iq — 1 



(6) 



where p is the spin- wave amplitude. Equation ([5]) also 
supports a bound state of <p(£) = ^sech£ for q = —i 



(uj = 0) [2fj,[21|. Equation ^ describes propagating spin 
waves without reflection, and takes an asymptotic form 
of tp (£ -> -oo) = pe iqi and tp (£ -> +oo) = -pj^e iq t. 
The spin wave maintains its amplitude and only captures 
an extra phase after passing through the DW. Interest- 
ingly, this phase shift is indeed observed in recent cal- 
culations [16H18| . The above result is very robust, and 
holds even with the extra Dzyaloshinskii-Moriya interac- 
tion [HE! Dm .(zx^) in Eq. ©. 

A very interesting consequence of the above results 
is schematically illustrated in the inset of Fig. 1: The 
magnons whose spins point to the left (opposite to the 
magnetization of the left domain) are injected into the 
DW from the left. The magnons transmit completely 
through the DW with their spins reversed (to the right). 
The change of magnon spins should be transferred to the 
DW, an all-magnonic STT. Thus the DW propagates to 
the left, opposite to the magnon propagation. One can 
also understand this result directly from Eq. ([T]). In the 



absence of damping, Eq. ([T]) can be cast as 

dm ^ „ d 

— — = m x Rm z z - —J, 
at oz 



(7) 
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also called spin-wave spin current 
nent of Eq. is conserved so that dtm z 
where J z is the z component of J. In terms of tp 



where J = Am x ^ is the magnetization current, 

The z compo- 
d z J z = 0, 
. Jz = 

(fd^V* ~ l P* a i l P) cos#o m the two domains. For the 
propagating spin wave (j6)), J z = —Ap 2 k in the far left 
of the wire (z — > —oo and 6q = 0), while J z = Ap 2 k in 
the far right (z — > oo and #o = n ), where k = q/A is the 
spin-wave vector in real space. The spin current changes 
its sign after passing through the DW, and results in an 
all-magnonic STT on the DW. Thus, in order to absorb 
this torque, the DW must propagate to the left with the 

2 

velocity V DW = —^-V g z, where V g = duj/dk = 2Ak is 
the group velocity. 

To test the validity of these findings in the realistic 
situation when both damping and transverse anisotropy 
are present, we solve Eq. (JlJ numerically in a one- 
dimensional magnetic nanowire. In the simulations, the 
time, length, and field amplitude are in the units of 
(7M S ) -1 , y/A/AI s , and M s , respectively, so that velocity 
is in the unit of r yy / AM s . If one uses the YIG parame- 
ters: M s = 0.194 x 10 6 A/m, K = 0.388 x 10 5 A/m, and 
A = 0.328 x lO" 10 A m 0, these units arc 1.46 x 10" 10 
s, 13 nm, 2.44 x 10 3 Oe, and 89 m/s. The wire length 
is chosen to be 1000 (from z = —500 to z = 500) with 
open boundary conditions and a transverse DW is ini- 
tially placed at the center of the wire. Spin waves are 
generated by applying an external sinusoidal magnetic 
field h (t) = ho sin (ftt) x of frequency and amplitude 
ho locally in the region of [—60,-55] in the left side of the 
wire. Thus, the spin wave (may not be monochromatic 
as explained later) propagates from the left to the right 
as illustrated in Fig. 1. We solve Eq. ([T]) numerically by 
using the standard method of lines. The space is divided 
into small meshes of size 0.05 and an adaptive time-step 
control is used for the time evolution of the magnetiza- 
tion. In terms of YIG parameters, the geometry of our 
nanowire is 0.65 nm x 0.65 nm in cross section and 13 
pm in length. The DW will move under the influence of 
the spin wave. The spatial-temporal dependence of m z 
is used to locate the DW center which, in turn, is used 
to extract the DW velocity. 

Below, we present our simulations for a set of realis- 
tic material parameters of YIG: a = 10 -5 and K± = 
2 x 10~ 3 . We present also the simulation results when 
both damping and transverse magnetic anisotropy are 
absent in order to show the quantitative effects of damp- 
ing and transverse anisotropy although the qualitative 
results are the same. Figure 2(a) is the numerical results 
of the spatial-temporal dependence of m z at ho = 1 and 
O = 0.75 (optimal frequency explained later) for YIG. 
The simulations show the following interesting results. 
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FIG. 2: (Color) (a) Density plot of m z in the z-t plane at 
frequency Q, = 0.75 for a = 10" 5 and K± = 2 x 10" 3 (YIG 
parameters). The center of the DW is initially at z — 0. The 
spin wave is generated in the region of z 6 [—60, —55]. (b) The 
frequency dependence of DW velocity. The red circles are for 
YIG parameters, and the blue squares are the results of the 
case without damping and transverse magnetic anisotropy. 



First, spin waves are generated by the external field. The 
spin waves propagate to both sides of the wire, resulting 
in the parallel strap pattern in the density plot of m z . 
Second, when the spin waves reach the DW at z = 0, 
the DW starts to move towards the left, opposite to the 
spin-wave propagating direction. The straight trajectory 
of the DW center before hitting the wave source indicates 
that the DW propagation speed is almost a constant. 
Third, the slopes of the spin-wave straps and DW tra- 
jectory tell us that the DW propagation speed is smaller 
than the spin-wave group velocity, a reasonable result 
that is consistent with our picture. It is also clear that 
there is no reflection when the spin waves pass through 
the DW in the presence of both damping and transverse 
magnetic anisotropy. This is highly nontrivial since it is 
not so clear from our earlier analysis. We will present 
further evidence for this finding. Figure 2(b) shows the 
frequency dependence of DW velocity at ho = 1 for both 
YIG parameters (circles) and the case without damp- 
ing and transverse magnetic anisotropy (squares). The 
error bars are smaller than the symbol sizes. The com- 
plicated and irregular frequency dependence of the DW 
velocity at low frequency is probably related to the ob- 
servation of the polychromatic spin-wave generation. At 
a large enough frequency (51 > 0.55), the excited spin 
wave is almost monochromatic with the same frequency 
as the oscillating field. These curves show that the DW 
propagation velocity is very sensitive to the microwave 
frequency. In fact, there exists an optimal frequency at 
which the DW velocity is maximal for a given set of pa- 
rameters. For the cases shown in the figure, the optimal 
frequencies are Q = 0.75 in the presence of damping and 
transverse magnetic anisotropy and a higher optimal fre- 
quency Q = 0.85 without them. 

The rcflcctionless property (total transmission) of the 
spin wave through DW can also be verified through quan- 
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FIG. 3: (Color online) Field dependence of spin-wave ampli- 
tude differences at two sites located in the opposite sides of 
the DW. Red circles are for the YIG parameters at the op- 
timal frequency (fl = 0.75), and blue squares are the results 
without damping and the transverse magnetic anisotropy also 
at its optimal frequency (fi = 0.85). 

titative analysis of spin- wave amplitude on the two oppo- 
site sides of the DW. If the spin wave is monochromatic 
(a sinusoidal wave) and passes through the DW without 
reflection, the difference of the spin-wave amplitudes on 
the two sides of the DW is around zero. We evaluate the 
spin- wave amplitude difference at z = —160 and z — 45 
at the same time, denoted as Sp. Figure 3 is the nu- 
merical results of Sp as a function of microwave field ho- 
Indeed, Sp is almost zero (green dashed line) both with 
(circles) and without (squares) damping and transverse 
magnetic anisotropy. Of course, for nonmonochromatic 
(sum of many sinusoidal waves) spin waves, the ampli- 
tudes on the two sides of the DW may be different at 
any particular time due to the complicated interference 
of waves with different frequencies. This is indeed the 
case for large ho, as shown by the oscillatory 5p around 
zero. The total transmission of spin waves through a DW 
is an important property because it results in a larger 
spin- wave spin current, and generates a larger magnonic 
STT. 

Figure 4(a) shows the ho dependence of the spin- wave 
amplitude p. It is almost linear at low fields both with 
(circles) and without (squares) damping and transverse 
magnetic anisotropy. The behavior is complicated at high 
fields, and a large error bar of p is observed, accompany- 
ing less regular and polychromatic spin- wave generation. 
This also results in a large fluctuation of p. The ho depen- 
dence of the DW velocity Vhw is shown in Fig. 4(b). It is 
nonmonotonic for the realistic situation with YIG param- 
eters (red circles) and almost quadratic for the case with- 
out damping and transverse magnetic anisotropy (blue 
squares). Although the field dependence of DW velocity 
is nonmonotonic, the relationship between the DW ve- 
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FIG. 4: (Color online) (a) Field dependence of spin-wave am- 
plitude for YIG parameters (red circles) at frequency Q, — 
0.75 and the case without damping and transverse magnetic 
anisotropy (blue squares) at frequency Q — 0.85. (b) Field 
dependence of the DW velocity for the same cases as those 
in (a). Inset: DW velocity vs the square of the spin- wave 
amplitude. Symbols are the simulation data, and solid lines 
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are ^-V s without any fitting parameters. 



magnonic STT could be made of magnetic insulators like 
YIG so that the Joule heating is, in principle, avoided. 
It is also known that the stray field is important in 
DW dynamics. Our results here are consistent with the 
OOMMF [II] simulations including this field. Remark- 
ably, these results are consistent with a phenomenologi- 
cal theory on thermomagnonic STT proposed by Kovalev 
and Tserkovnyak [29j | . 

In conclusion, we proposed an all-magnonic spin- 
transfer torque mechanism for magnetic domain wall ma- 
nipulation in nanowircs. This spin-transfer torque can 
effectively drive a DW to propagate along the wires. 
The propagation speed is sensitive to both microwave 
frequency and its amplitude. There is an optimal fre- 
quency, order of the usual ferromagnetic resonance fre- 
quency, at which DW propagating speed is the fastest. 
All-magnonic STT should have advantages over its elec- 
tronic counterpart on energy consumption as well as on 
the spin-transfer efficiency. It also opens the door for 
using magnetic insulators in spintronic devices. 

This work is supported by Hong Kong RGC 
Grants (No. 604109, No. RPC11SC05, and No. 
HKUST17/CRF/08). 



locity and the spin-wave amplitude is much simpler. As 
shown in the inset of Fig. 4(b), the DW velocity Vdw 
is almost quadratic in p both with (circles) and without 
(squares) damping and transverse magnetic anisotropy. 
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We also plot ^V g (solid lines) without any fitting pa- 
rameters, where V g = 1.48 at £1 = 0.75 for the YIG case 
and V g = 1.61 at £1 = 0.85 in the absence of damping and 
transverse anisotropy, and p is calculated numerically. 
Although there is no reason why the early velocity for- 
mula derived under the approximation of zero damping 
for uniaxial wire and small spin-wave amplitude should 
be applicable to the realistic case when both damping 
and transverse magnetic anisotropy are presented, the 
theoretical formula is, in fact, not too far from the nu- 
merical data for both cases. Of course, it should not be 
surprising for the deviation at large p since quadratic p 
dependence of Vdw is derived based on the conservation 
of the z component of angular momentum that does not 
hold for the generic cases. Also the main purpose of the 
current study is to demonstrate the principles rather than 
the exact mathematical expression of DW velocity which 
can be the subject of future studies. 

Most studies [T3-E1 HI US H3] of magnonic effects 
in nanomagnctism so far are about the conversion of 
magnon spins with electron spins. Very often, it goes 
through the Seebeck effect that involves both thermal 
and electronic transport. Thus, like usual electronic 
STT, devices based on these effects must also contain 
metallic parts so that Joule heating shall be present. 
In contrast, the magnonic STT presented here does not 
require electron transport. Devices based on this all- 
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